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1. HQET: Physical Motivation 

Heavy quark effective theory (HQET) is the Umit of QCD where the heavy 
quark mass, mg, goes to infinity with its four- velocity v^^ held fixed. This 
limit is appropriate for physical systems where there is a heavy quark inter- 
acting with light degrees of freedom that typically carry four-momentum 
much less than the heavy quark mass. Consider, for example, a meson with 
Qq flavor quantum numbers. (Here Q denotes a heavy quark and q a light 
quark.) The size of such a meson is set by the nonperturbative scale of the 
strong interactions, r 1/Aqcd- Hence, by the uncertainty principle, the 
typical momentum carried by the light degrees of freedom is, pg ^ Kqcd- 
If an amount of momentum 5pi is transferred to the heavy quark the change 
in its four- velocity [pq = rngv) is 6v ~ Spi/mQ — > 0, in the infinite mass 
limit. That is why, for the study of such hadrons, an appropriate limit of 
QCD involves taking ttiq oo and holding the heavy quark's four-velocity 
fixed. As mg oo the strong interactions of the heavy quark become in- 
dependent of its mass and spin. Consequently there are new spin-flavor 
symmetries that arise in this limit. These approximate symmetries endow 
HQET with predictive power in the nonperturbative regime. 



2. HQET: The Effective Field Theory 

The part of QCD Lagrangian containing a heavy quark is 



where Q is the heavy quark field and D denotes a covariant derivative. We 
want to take the limit mg — > oo with fixed four- velocity v^^. To do this 
one scales out the rapidy varying part of the heavy quark field's space-time 
dependence by writing. 



jCqcd = Q{ip - mQ)Q. 



(2.1) 



(2.2) 
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Putting this into the QCD Lagrangian 

jO.QCD QvirriQii + ip- mQ)Qy = QytpQy 



up 



iv ■ D + if 



1 



(2.3) 



(2.4) 



2 y V 2 

So in this hmit the QCD Lagrange density becomes [1] 

The HQET Feynman rules follow from the above Lagrangian. The heavy 
quark propagator is 

^ (^) ■ <-) 

and the heavy quark gluon interaction vertex is 

-igvf'T^, (2.6) 

where is a color generator. This can also be derived by taking the 
appropriate limit of the QCD Feynman rules. Writing the heavy quark 
momentum as, pQ = rnqv + k , the limit is taken by neglecting the residual 
momentum k in comparison with rnqv. For the propagator this gives, 



imQ{i; + 1) 
2'mQV ■ k + is 



V ■ k + ie 



(2.7) 



For the the heavy quark gluon vertex. 



7m 



-igT\^. 



(2.5 



Here we used the fact that the vertex always appears between propagators 
so factors of ij) + l)/2 can be inserted. 

HQET is a theory of the heavy quark. The field Qy destroys a heavy 
quark but it docs not create the corresponding anti-quark. There is no pair 
creation in HQET. This follows straightforwardly from the Feynman rules. 
Notice that the propagator only has one pole in So a closed loop of 
heavy quark propagators will have all the poles in the zero component of 
the loop momentum integration below the real axis. Thus the integration 
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vanishes. (When the contour of this integration is completed in the upper 
half plane it encloses no poles.) 

We have motivated physically the fixed heavy quark velocity superse- 
lection rule. Lets now try to motivate it mathematically. To derive the 
effective field theory for a heavy quark of four-velocity vi and a heavy 
quark of four-velocity V2 write, 

Q ^ g-m,.i.XQ^^ _^ ^..rUQV^-.Q^^^ (2.9) 

where ^iQvi = Qvi and i)2Qv2 = Qv2- Putting this in the QCD Lagrangian 

gives, 

i^QCD + £j,2 -t- (2-10) 

The cross terms that mix the heavy quark fields of different velocity is 

Cv,v2 = e-'^Q^^^-^^^^Q^^lmg^a mQ)Qy, + h.c. 

= e-^™o(''i-''=)^Q„,ipQ„, + h.c. (2.11) 

This contributes a negligible amount to the action in the large mass limit 
because of the rapid variation of the exponential. Similarly the measure in 
the path integral factorizes 

[dQ] ^ [dQ,,][dQ,,], (2.12) 

since the two types of modes are orthogonal. So we see that 

Chqet = Qviv ■ DQy. (2.13) 

V 

For a given four-velocity v the HQET Lagrangian is independent of 
the heavy quark's mass and spin. However, it does depend on the heavy 
quark's velocity. All heavy quarks with the same four-velocity interact in 
the same way. If there are N heavy quarks at velocity v then the theory 
has a SU{2N) spin- flavor symmetry [2]. Note that this symmetry is a little 
unusual since it relates states at the same velocity but different momentum. 

Suppose we want to describe an anti-quark instead of a quark. We start 
at the same place except write 

Q = e+imQV.XQ_ ^Q_^_Q_^ (2.14) 

and we anticipate Qy will be the field that creates an anti-quark. Putting 
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this into the QCD Lagrange density 
>C = Qv{-mQ'^ - mq + iIl))Qv 

= -Qyiv ■ DQy. (2.15) 

Heavy anti-quark Fcynman rules follow from this or from taking the ap- 
propriate limit of the QCD Feynman rules. Writing the anti-quark four- 
momentum as, pQ = —rnqv — k, the heavy anti-quark propagator is 



1-^ 



(2.16) 



V ■ k + ie \ 2 

while the heavy anti-quark vertex is 

+ig{T^fv^. (2.17) 

Thus the heavy anti-quark Feynman rules are identical to those for the 
heavy quark except for the replacement {—T^Y". This is expected 

since the anti-quark is a color 3 and the color generators acting on this 
representation are {—T^)'^ = . 

3. 1/toq Corrections to the HQET Lagrangian 

There are l/mq corrections to the HQET Lagrangian. To find these lets 
reexamine the derivation of HQET. For this we wrote, 

Q = e-^^Q^-^Q, ^Q„ = Q„. (3.1) 

Factoring out the phase factor is not an approximation it is just a field 
redefinition. The approximation comes in putting ^Qy = Q„. So to go 
beyond leading order in l/mq we make the most general decomposition 

■j/iQy = Qy, fxv = -Xv (3.2) 



Now the decomposition is completely general. The field x corresponds 
partly to the anti-quark degrees of freedom. Putting this into the QCD 
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Lagrange density gives 
C = Q{ip - mQ)Q 

= [Qv + xv]{mQ{i) - 1) + + xv] 

= Qyiv -pQv - Xv{iv - D + 2mQ)xv 

+ QvipxXv+Xvi-p±Qv (3.3) 

It is convenient to introduce a notation for the part of a four-vector per- 
pendicular to the heavy quark's four- velocity v. For any four- vector V its 
perpendicular part V± is defined by 

Vj;; = V^-{vV)v^', (3.4) 

and V ■ V± = 0. Because Qvl^Xv = XvfQv = we have replaced p in the 
last two terms by^j,. 

The field Xv has a mass 2mQ and consequently we can integrate it out 
of the theory. At tree level this corresponds to solving the equations of 
motion to express it in terms of Qy. The field equation 

{iv - D + 2mQ)Xv = ipi.Qv 
implies that 

= (■ nio (3.5) 

(iv-D -V 2mQ) 

Putting this back into the QCD Lagrange density gives 

£ =Qy{ivD + ip^ (2mQlivD) ^-l} Qv 

= Q^iv ■ DQ, + ^^Q.ipx^jp^Qy + O (^) . (3.6) 

The relationship between the heavy quark field in QCD and in the effective 
theory is 



Q 



-rriQV 



2mQ 



Qv + 0\^]. (3.7) 



You will often see these formulas written without the _L subscript. That is 

because if one works at order 1/toq the equation of motion v ■ DQ^ = 
allows one to write D± as a full D. Expanding in powers of l/niQ, 

£ = £o+Ci + ..., 

where the leading term is the HQET Lagrange density 

Co = I^HQET = Qviv ■ DQy. (3.8) 
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To simplify Ci write, 

p±p± = rrD^f.D±, = i{r,r}D±^D^^ + ^[r,r]D^D^ 

= Dl + ^[r,r][D^,D.]- (3.9) 

Since this is sandwiched between Qy and Qv the Z)_l's can be replaced 
by D's in the anticommutator term. Using cr^'^ = ^[y^,^"^] and igG^i, = 
[Df^,D,^], the above becomes, 

P±P± = Dl + ^a^-G^,. (3.10) 

Using this in eq. (3.6) gives the order l/mg term in the Lagrange den- 
sity [3]. 

- n ni^^ 

jCi = -QvTT^Qv - gQv-^ Qv (3.11) 

The first term in £1 is the heavy quark kinetic energy. It breaks the flavor 
symmetry but not the spin symmetry. The second term is a magnetic 
moment interaction /Xq - Be. It breaks both the spin and flavor symmetries. 

This is the Lagrange density at tree level. The coefficients change when 
one includes perturbative connections. Most importantly the operator 
Qv'^fivQ^'^Qv requires rcnormalization so its coefficient develops a factor 
a(/i), where the subtraction point dependence cancels that of the operator. 
The tree level matching of QCD onto HQET that we have performed so 
far implies that 

aimq) = 1 + 0(qs(toq)). (3.12) 

In dimensional regularization with minimal subtraction the kinetic energy 
operator does not get renormalized to all orders in perturbation theory 
and its coefficient stays equal to unity. One way to understand this is as a 
consequence of reparametrization invariance. 



4. Reparametrization Invariance 

A heavy quark's four-momentum can be written as the sum of rngv and 
the residual momentum k, 

PQ = rnqv + k. (4.1) 

This decomposition is not unique. Typically k ^ 0{Kqcd) and a small 
change in v of order I^qcd I'mq can be compensated by a change in k of 
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order Aqcd- Explicitly 
V ^ V + e/iriQ, 

k^k-s, (4.2) 

leaves the heavy quark's four-momentum unchanged. Since the four- 
velocity satisfies = 1 we must, at linear order in e, choose 

v-s = 0. (4.3) 
We also want to preserve the constraint 

^Qv = Qv (4.4) 
Therefore if 

Qv Qv + SQv, (4.5) 
then SQv satisfies 

^ + ^]{Qv+SQ^) = Q^ + SQ^, (4.6) 

which implies that 

(1 - i!)SQv = -^Qv (4.7) 

Hence we can choose 5Qy = 2^Qv' which satisfies -^SQy = —SQ^. There- 
fore the freedom to choose different heavy quark velocities v implies that 
the theory should is invariant under [4] 

V ^ V + e/ruQ 

Under this transformation the leading order Lagrange density 

jCo = Qviv ■ DQy, (4.9) 
transforms to 

implying that 

SjCo = —Qv{ie ■ Di_)Qy. (4.11) 
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Here wc used Q^^Qv = which follows from the constraint e ■ v = 0. The 
change in £i is 

SCi ^ -Qvi2^-^Qy = — ^—Q^ie ■ D^Q^. (4.12) 

Therefore jCq + CL\ is invariant under the transformantion in eq. (4.8) 

provided the coefficient of the kinetic energy term is 1 . Hence if a regulator 
is used that preserves invariance under reparametrization transformations 
the kinetic energy cannot get renormalized. 

We have been developing formalism for a h^qcD l^^-Q expansion of phys- 
ical quantities. The quarks that are heavy enough for this to be useful 
are the charm, bottom and top which have masses mc ^ 1.4 GeV, 
4.8 GeV, and mt ~ 175 GeV. However rat is so heavy that it decays before 
forming a hadron {t ^ h + W, with a width of a few GeV). So the tools I 
am developing here are only useful for the charm and bottom quarks. 



5. Spectrum of Hadrons Containing a Single Heavy Quark 

In the mq — > oo limit hadrons containing a single heavy quark are clas- 
sified (in their rest frame) not only by their total spin (i.e., total angular 
momenta) S but also by the total spin of the light degrees of freedom [5], 

S = S, + Sq. (5.1) 

For any angular momenta J its square is = j(j + 1). Since sq = 
\ hadrons with a given spin of the light degrees of freedom come in 
degenerate doublets with spins 

s = st±^. (5.2) 

In the case of mesons with Qq flavor quantum numbers two doublets have 

been observed in the Q = c and Q ~ b cases. The ground state doublet has 
S£ = ^ and negative parity. For Q = c these spin and spin 1 mesons are 
the D and D* while for Q = b they are the B and B* mesons. The other 
observed doublet has = | and positive parity. For Q = c the spin 1 and 
spin 2 members of this doublet are called the Di and D2 respectively. 

There is a very useful phenomenological model that describes with sur- 
prising accuracy the hadronic spectrum. It is called the nonrclativistic 
constituent quark model. In this model the light u, d quarks have masses 
of 350 MeV and the strange quark has a mass of 500 MeV. The quarks are 
nonrelativistic and interact via a potential that is linear at large distances. 
The quarks in the nonrelativistic quark model are not the ones destroyed 
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and created by the fields in the QCD Lagrangian. They are quasi parti- 
cles and their large masses arise from nonperturbative strong interaction 
dynamics. In this model the ground state Qq meson has the light anti- 
quark in an S-wavc and the total angular momentum of the light degrees 
of freedom comes from the light quark spin. Consequently this doublet has 
si = \ and negative parity coming from the intrinsic parity of an anti- 
quark. The first excitations correspond to giving the light anti-quark one 
unit of orbital angular momentum L = 1. Combining this orbital angular 
momentum of the light degrees with the spin of the anti-quark gives the 
possibilities. ,S£ = ^ and |. The parity is now positive because of the unit 
of orbital angular momentum. As was mentioned before only the | doublet 
has been observed. 

Heavy quark symmetry and a little phcnomenological "lore" helps pro- 
vide us with an understanding of why the other, si = ^, doublet is not 
observed. The excited mesons in the Si — ^ and | doublets decay to the 
ground state mesons with emission of a single pion. The available phase 
space is about 450 MeV and decay to more pions is probably suppressed. 
Parity invariance implies that the one pion decay occurs in an even partial 
wave, either L = or 2. The decays D2 Dn and D2 D*tt, are L = 2 
while decay Di D*n could occur in either L = 2 or L ~ 0. But in the 
rric — > 00 limit the £>i and total widths must be equal, so the L = 
amplitude is forbidden. On the other hand the S£ = ^ doublet of excited 
charmed mesons (denoted and Djf) decay to D^*^Tr in an S'-wave. Am- 
plitudes that go through high partial waves are smaller than those that go 
through lower partial waves (they are suppressed by (pTr/lGeV)^^"'"^) and 
so the Di and D2 are narrow having widths around 20 MeV. The members 
of the excited se = ^ doublet decay in an S'-wave and are expected to be 
too difficult to observe because they are quite broad (i.e., widths greater 
than ~ 100 MeV). 



6. V- A = Gauge 

Calculations in HQET can be performed in almost any gauge. However 
in V ■ A = gauge HQET perturbation theory is singular. Consider tree 
level Qq elastic scattering in the rest frame v = v,. = (1,0). In HQET 
an on-shell heavy quark has a four-velocity v and a residual momentum k 
that satisfies v ■ k = 0. Suppose the initial heavy quark has zero residual 
momentum and the final has residual momentum k = (0, k). The tree level 
Feynman diagram shown below in Fig. (1) gives the heavy quark light 
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quark scattering amplitude 



M = g'^u{vr)T^u{vr):^u{k)T'^MO) (6.1) 

in Feynman or Landau gauge (here u{k)^u{0) = was used). However in 
V ■ A = gauge the gluon propagator is 



fc2 



■ te 



V ■ k [v ■ ky 



(6.2) 



The heavy quark kinetic energy cannot be treated as a perturbation in this 
gauge because this imphes that v-k = Q and the square brackets are ill de- 
fined. Including the heavy quark kinetic energy in the Lagrangian the resid- 
ual momentum of the outgoing heavy quark becomes = (k^/2mQ,k) 
and V ■ k = k^/2TOQ is not zero. At leading order it is the QQA vertex 
from the heavy quark kinetic energy that contributes to elastic Qq scatter- 
ing in V ■ A = gauge [6]. The Feynman rule for the Qv{k') Qv{k)+ 
gluon vertex arising from an insertion of the kinetic energy operator is, 
{g/2mQ){k_L + k'^_)^ = {g/2mQ){k + k')^ + {g/2mQ)v ■ [k + k')v^. The part 
proportional to t;^ doesn't contribute, and since u{k)^u{0) = only the 
second term in the numerator of the propagator matters. For large tuq it 
reproduces the amplitude M. above. 



Fig. 1. Qq scattering via one gluon exchange. 

For the calculation of on-shell amplitudes in t; . A = gauge the heavy 

quark kinetic energy must be considered as leading and it is the QQA vertex 
from this operator that gives rise to the leading Qq scattering amplitude. 
However, for off-shell Green functions this gauge can be used in HQET, 
with the kinetic energy treated as a perturbation. Since anomalous di- 
mensions follow from such Green functions we can conclude, for example. 
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that the anomalous dimension of the operator Qv^Qv vanishes, without 
performing any calculation. There is another way to understand this. The 
operator QyVQ^ is a charge density of the heavy quark spin-flavor symme- 
try and therefore cannot be renormalized. 



7. Mass Formulae 



The masses of hadrons containing a single heavy quark can be expanded 
in powers of l/ruQ. In the rest frame v = Vr = (1,0) 

-Mjj(Q) = mg + A H 



2 • 2m 



Q 



+ a[fi) 2.4^Q +•••' ^^■^) 

where the ellipses denote terms suppressed by more powers of l/mq. A 
is the energy of the light degrees of freedom in the mg oo limit. It is 
the same for both members of a spin symmetry doublet but of course is 
different for different doublets. The kinetic energy term is also the same 
for both members of a given doublet and it is conventional to write 

2X, = -{H(,Q)\Q,^DlQ,JH(Q)), (7.2) 

Ai is independent of the heavy quark mass and is of order Aq(jjj. The 
matrix elements in cq. (7.1) arc taken in HQET with the hadron states 
having zero residual momentum. There is a factor of two in the denomina- 
tor because of the HQET normalization convention, {H'^9\k')\Hi'^\k)) = 
2v° {27r fS^.y,S^{k - k'). Since Qv,(JafiG°''^Qy^ = Q.^a ■ B^Q^v tlio matrix 
element of this operator is proportional to Sg • S^. It is conventional to 
write 

(Sg • S,)16A2 = -a(M)««)|g..5^Ta/3G"'^Q.J^^i?)). (7.3) 

We can simplify this using 

Sg • = (S^ - Sg - S|)/2 

= (s(s+l)-3/4-s^(s^ + l))/2. (7.4) 
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For the higher spin (s+) member of a doublet se = s+ — 1/2 and the above 
is 

2(Sq • S^)+ = s+{s+ + 1) - 3/4 - (s+ - l/2)(s+ + 1/2) 
= sl+s+- 3/4 - 4 + 1/4 

= i(s+-l) = i(2s_ + l), (7.5) 

while for the lower spin (s_) member of a doublet 

2(Sq • S^)_ = s_(s_ + 1) - 3/4 - (s_ + l/2)(s_ + 3/2) 
= sl+s-- 3/4 -si- 2s- - 3/4 

= -^-(2s_+3) = -l(2s+ + l). (7.6) 

So 

(SQ.S,)±=±in^. (7.7) 

where nip = (2s:p + 1) is the number of spin components. The matrix 
element A2 is defined 

A2 = T^a(/z)(i/;g|g.,ffaa0G"'^g.Jij{g), (7.8) 

and the mass formula is 

m^(Q) = mg + A - ± + . . . . (7.9) 

^(±) 2mQ 2mQ 

Note that the spin averaged mass 

^H(Q) = ^, (7.10) 

n+ + n_ 

is independent of A2. Ai is independent of the heavy quark mass and A2 
has a weak logarithmic dependence on mg. The difference in mass between 
the two members of a doublet determines A2. 

{n-+n+) ^ {2si + 1) , 
rn„(Q) - m„(Q) = — A2 = A2. (' -ll) 

"{+) "(-) 2mQ ruQ 

For the ground state multiplct of mesons with se = ^ the B* — B splitting 
implies that A2(mb) = 0.12 GeV'^ while the D* — D splitting gives Ac(toc) = 
0.10 GeV^. Here (and throughout these lectures) I use rub = 4.8 GeV and 
rric = 1.4 GeV. The charmed states in the sj'^ = I"*" multiplct arc the 
1)2(2460) and Z?i(2420). The mass splitting between these states implies 
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that for this multiplet AJ, = 0.013 GcV^. In these lectures I use the notation 
A, Ai,A2, mp for the ground state meson multiplet and A' , A'^ , A'^ , mp for 
the excited s^^ = |^ doublet. 

The values of A, Ai arc not well known. A determination from the lepton 
energy spectrum in inclusive semileptonic B decay gives A ~ 0.4 GeV, Ai ~ 
—0.2 GeV^, but the uncertainties are very large [7]. 

Using the mass formula above it is straightforward to deduce that [8] 

-^,_-^^m,{m'B-i^B)-mM'n) ^^ ^^ GeV 

a; - Ai = ^^cm,[{rh'^-mB)-{m'^-fnn)] ^ 
mb - rric 

Here I used =5.73 GeV [9]. There is considerable uncertainty in this 
value because the peak of the 3*71 mass distribution measured at LEP may 
not correspond to the narrow |^ doublet. The above formula for A' — A 
holds up to corrections of order ^qcd/^q since the dependence on Ai 
cancels out between the two terms in the numerator. The extraction of 
A' — A is not very sensitive to the numerical choices for nib^c- 



8. Covariant Representation of Fields 

The fields in the ground state meson doublet are denoted by Pi'^^ and 
P*i'^^ . The vector meson field satisfies v'^P^Jf^^ = 0. It is convenient to 
combine the P^^^ and P^^j^^ fields into a single object that transforms 
in a simple way under heavy quark spin transformations. The mesons are 
composed of a heavy quark Q and light degrees of freedom as is the bispinor 
So the object we want should transform in this way under Lorentz 
transformations. The obvious choice is [10,11], 

^(Q)^(l+^[^^p,,(Q)_^3p(Q)]. (8.1) 

The bispinor field satisfies 

i,Hi'^'> = iJ^Q) and = ~Hi'^\ (8.2) 

and under heavy quark spin transformations 

Hi^^ ^ D{R)QHi'^\ (8.3) 

The matrix D{R)q is the usual Dirac four-component spinor representation 
of rotations and it satisfies ^^D\R)q^^ = D{R)q^. The 75 in eq. (8.1) is 
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dictated by parity, under which Hy'^^ transforms as 
where 

xp = -x), VP = {v°, -v). (8.4) 
Note that 

HP = [rP:P' + l.P^'''']^^ (8.5) 
and that 

TrHi^^^Hi^^ = -2Pi«)t^(Q) + 2P;^«)+P;W)'^, (8.6) 
is invariant under Lorentz transformations and heavy quark spin symmetry. 

9. Semileptonic B Z)(*'ei>e Decay in the mq oo Limit 

In semileptonic K TTei>e decay it is conventional to use the variable = 
{Pe+PuJ"^ to describe the kinematics of the decay. For B D^*^eve decay 
= iPe+Pvef' = (PB -PdmY is in the range < < (ms - rrijyM)'^. 
The maximum value of corresponding to the zero recoil kinematic point 
where the D^*^ is at rest in the rest frame of the B meson. The minimum 
value of q^ occurs when the D* is recoiling maximally. Even though q^ can 
be of order ^ the momentum transfer to the light degrees of freedom is 
much less than the heavy quark masses. The momentum transfer to the 
light degrees of freedom is of order 

Qught ~ i^QCDv' - AQCDvf = 2A^cfl(« • ^' - !)• (9-1) 

For B D'^*'>eve decay a better variable than to use is the dot product 
of B and four-velocities w = v ■ v'. It is related to q"^ by, 

, ml + m^(.) - 
w = v-v = . (9.2) 

2mBmD(.) 

For these decays 1 < w ^1.6 and hence qfig^t is small compared to the 
heavy quark masses over the entire phase space. Consequently HQET is the 
appropriate starting point. The invariant matrix element for the semilep- 
tonic decay is (since the electron and neutrino don't interact strongly) 

M = ^V,t{D^*\p')\n^{l - j,)b\B{p))[u{pe)r{l - J5)v{p.j]. (9.3) 
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All the complicated nonpcrturbative strong interaction physics is in matrix 
elements of the vector current = cj^b and axial vector current = 
c^iil^b. These matrix elements can be written in terms of Lorentz scalar 
form factors and it is convenient for comparisons with the predictions of 
HQET to write these form factors as functions oi w = v ■ v' , where the 
B meson has four-momentum p = uibV and the D^*^ meson has four- 
momentum p' = mjj(,)v'. When wc match onto HQET wo use this choice 
of four-velocities for the heavy quark fields and the heavy meson states. 
This corresponds to setting the residual momentum of the B and D^*^ 
states to zero. In terms of Lorentz scalar form factors the matrix elements 
are 

{D{p')\V>^\B{p)) 



{D*{p' e)\V^^\B{p)) 

{D*{p' e)\A^\B{p)} 



= h+{w){v + v')^ + h-{w){v - vY 



-hA3{w)is* ■v)v^. (9.4) 

There are several features of these formulas that need explanation. 
Firstly note that for the B D case there is no axial current matrix 
element. This is a consequence of parity invariance of the strong interac- 
tions. Suppose we write 

{D{p')\A"\B{p)) ^ ^^^^^^^ ^ ^ hA_{w){v - v'Y, (9.5) 

under parity 
P-^AP = A, 

P-^A°P = -A°, (9.6) 
and 

P\Dip)} = -\D{pp)} (9.7) 
where 

PP = {P°,-P)- (9.8) 
Now 

{D{p'p)\A'^\B{pp)) = hA4w){vp+v'pr + hAAw)iyp-v'pr 

= {PD{p')\A>^\PB{p)) = {D{p')\p-^A'^P\B{p)). (9.9) 
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This implies that 

- [hA+ {w){v + v'f + hA_ {w){v - v'f\ 

= hA+{w){v + v')° + hA_{w){v - v'f, (9.10) 
and 

[hA+ iw){v + v'Y + hA_ {w){v - v'Y] 

= -hA+iw){v + v'y -hA_{w){v-v'y. (9.11) 
Consequently 

hA^ = hA_ = 0. (9.12) 

Similarly the i in the matrix element of the vector current between B and 
D* is from time reversal invariance of the strong interactions. Finally the 

factor of ^Jvn.BfTi.jji,) in the denominator occurs because of our normaliza- 
tion of states. In full QCD wc use the standard normalization 

{H{p')\H{p)) = 2E{2T:f5\p - p'), (9.13) 

which has a factor of mn from E. With y^niBmjjM in the denominator 
the matrix elements on the left hand side become independent of mg in 
the large mass limit. 

The first step in using heavy quark symmetry to constrain the form 
of the matrix elements is to match the weak current onto an operator in 
HQET. Neglecting as{mc,b) and ^QCD/fnc,b corrections this matching is 

C7p75& = Cy>jfj,j5by. (9.14) 

Operators Cy'Tby transform in a particular way under heavy quark spin 
symmetry. If we pretend that T — > Dc{R)TD~^{R)b then the operator is 
invariant. For B — ^ matrix elements we want to represent the operator 
Cv'^by in terms of B and Z?*^*-* fields constructed so that it transforms in 
the same way as the quark operator. This gives 

c„.r6„ = TrHl,';^TH^''^X, (9.15) 

where Xap is a general matrix written in terms o{v,v' the gamma matrices 
and the identity matrix, 

X = Xo + Xiij + X2f +X30. (9.16) 

But since i>Hi''^ = H^^'^ and iJ^fV = each term in eq. (9.16) has 

the same effect. We can equivalently write, 

X = -$(«;), (9.17) 
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where £,{w) is usually called the Isgur-Wisc function. Taking the B D^*'> 
matrix element of eq. (9.15) implies the following relations between the 
form factors [2] 

h+{w) = hv{w) = hAiiw) = liA^iw) =■ i{w) 



At zero recoil c„r6b is a generator of spin flavor symmetry and its matrix 
element is known. This implies that [2,12] 



The HQET Feynman rules are close to the full QCD Feynman rules 

for residual momentum small compared with mq. However, in loop inte- 
grations the residual momentum gets arbitrarily large and for regions of 
integration where k > mg the HQET results are very different from those 
in full QCD. Fortunately QCD is asymptotically free and therefore these 
differences can be taken into account by adding perturbative corrections 
to the coefficients of HQET operators. The perturbative corrections to the 
matching of the weak currents onto HQET operators do not cause a loss 
of predictive power since dimensional analysis dictates that the HQET op- 
erators that occur in the matching condition are of the form c-v'Tby. The 
B D'^*^ matrix elements of any operator of this type (not just those oc- 
curing in eq. (9.14) where T is 7^ or 7^75 ) are determined by ^. I will not 
discuss such perturbative corrections in these lectures. An excellent review 
of this subject, with references to the original literature, can by found in 



10. Luke's Theorem for the l/ruQ Corrections 

There are two sources of l/mg corrections to the form factors, corrections 
to the currents and corrections to the states. Neither of them change the 
matrix elements of the weak currents at zero recoil. This result is important 
for the determination of the 6 ^ c element of the Cabibbo-Kobayashi- 
Maskawa matrix from exclusive B decays and is called Lukes theorem [14]. 
The corrections to the currents arise from the l/mq term in 



h-{w) = hA-iiw) = 0. 



(9.18) 



^(1) = 1- 



(9.19) 



[13]. 



Q = e 



—imQV-x 



1 + 



Q 



(10.1) 



2mQ 



and those from corrections to the states arise from £1 in eq. (3.11). I will 
discuss both types. For the currents including the order 1 / niQ term in eq. 
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(10.1) gives 

cTb = a,,Th, + c„, f-^r + 6,, (10.2) 

\ znic 2mb I 

for the relationship between the currents in full QCD and HQET. For the 

B — > D^*^ matrix elements of the 1/toq terms 

Cyn'D^TK = TrH[,';'>THi''^Ml''^ 

c^'ViD^K = TrHl;^TH^^'>M^!'\ (10.3) 
where the general form for the matrices M^i*^' is 

(^; + «')^ + d''H^-^')^-d'^^. (10.4) 

Here and ^g*^^ are functions of w. But we have the following 

identity valid for a i? ^ D*^*-* matrix element 

iD^{c^,TK) = K{v' - v)^,i{w)TrHi'}THi''\ (10.5) 

which implies the relations 

^f+eW=0. (10.6) 

Eqs. (10.6) allow us to express ^j'^ in terms of The leading order 
equations of motion imply that 

v'^C^'iD^'Tby = 0, and Vf,Cy,TiDf,by = 0. (10.7) 

This gives the following relations between the form factors 

+w)+ ^i''^ = (10.8) 

il + w)+ d*"^ il-w)+ = 0. (10.9) 

So all the current corrections can be expressed in terms of A, £,{w) and 
£,3{w). At zero recoil ^_ has no effect on the matrix element. At this 
kinematic point the above relations become 

2d'^(l)+d'^(l) = 0, 

2d'')(l)+?f (1)=0. (10.10) 
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Eqs. (10.10) imply that at zero recoil (where v' = v) Mj^ ' is proportional 
to 7^ + Vn which vanishes in the trace of eqs. (10.3). Hence we conclude 
that the order l/mg corrections to the current do not contribute to the 
B —> D^*^ matrix elements of the weak currents at zero recoil. 

Now we consider the Lagrangian corrections. These are corrections to 
the relation between heavy meson states in full QCD to those in HQET. 
They are represented by the HQET matrix element of 

i J d'^xT{Ci{x)icy>rK){0)}. (10.11) 

There are two types of terms in Ci. One is the heavy quark kinetic energy. 
It does not break the spin symmetry and so its effects can be absorbed into 

a redefinition of the Isgur-Wise function, £,{w). Next consider the charm 
quark chromomagnetic term. For B — s- D^*^ matrix elements, 

i [ d^xa{fi)T{{cy,ga„0G"f^Cy>){x){cy>Tby){O)} 



where 

X"'^ = iX2iw){v"^^ - A") + X3{w)a"^. (10.12) 
Since 

(^i J d''xa{^i)T{{b,gaa0G''^'h){x){h^c,,){O)}^^ 

= -iJ d*xain)T{{bygaa^G'"^by){x)ic,>rb^)iO)}, (10.13) 

the bottom quark chromomagnetic correction is also characterized by X2 
and X3. Only X3 can contribute at zero recoil. 

These Lagrangian corrections also preserve the normalization of matrix 
elements at zero recoil. The 1 /mq corrections change a heavy meson state 
to |p(Q)(*)) + (e/mQ)|5'('3)(*)) where the new state |S'(<3)(*)) is orthogonal 
to the niQ — > 00 state |pW)(*)), and e is a quantity of order Aqcd- At 
zero recoil c^Tby is a charge of the heavy quark effective theory. It takes 
\p(b){*)) into |PW(*)) and consequently (5(^)(*)|c„r6„|P(*')(*)) = 0. So 
there are no AgcD/mQ corrections at zero recoil. In other words even 
after absorbing kinetic energy correction into ^ we still have normalization 
condition ^(1) = 1 and furthermore ^3(1) = 0. 
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11. Kinematics of Inclusive B Decay 

Semileptonic B-meson decays to final states involving a charm quark arise 
from matrix elements of the weak Hamiltonian density 

Hw = ^VcbCji'il - 75)&e74l - ^5)1^^- (11-1) 

In 3-body exclusive semileptonic decay one looks at fixed final states, like 
Devg. The differential decay distribution has two independent kinematic 
variables which can be taken to be Eg and E,^^ . The decay distribution 
depends implicitly on the mass of the final hadronic state which is treated 
as a constant. In inclusive decays one ignores all details about the final 
hadronic state Xc and sums over final states containing a charm quark. In 
addition to the usual kinematic variable i^e and E,j^ that occur in exclusive 
semileptonic decay there is a third variable which will be chosen to be q^, 
the invariant mass of the virtual W^-boson. 

The differential decay distribution for inclusive semileptonic decay is 

■5{E,-pI)5{E,^-pI^) 

•E E ^^^^|^^(2-r^^bB-(Pe+P.J-pxJ, (11.2) 

\ spina J 

where we have used formula, cPp/2p'^ = d'^pS{p^ — m^). We work in the 
rest frame of the B meson. After summing over final hadronic states Xc 
the only relevant angle is that between electron and neutrino. Integrating 
over the direction of neutrino gives An Therefore 

d^Ped^P., = i'ln)i2TT)\pc\^d\pe\\p.fd\p,Jdcos0. (11.3) 

Using 6{Ee—Pe) and 6{Eu^ —plj to do integrations overpg andp°, 6{pl) = 
5{El — |peP) to do integration over \pe\, S{plJ = 5{El^ — |pi/J^) to do 
integration over |p^J and 5{q^ — {pe + PuJ^) = 5{q^ — 2E^E^^{1 — cos6)) 
to perform integration over cos 9 we have 

dT V- \{Xceu,\Hw\B)W ^^ r ^ ^ ^ 



dq^dE^dE^^ ^ f;^ Stub 

\ spins ) 



(11.4) 
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The square of the weak matrix element can be factored into a leptonic 
matrix element and a hadronic matrix element, since leptons do not have 
strong interactions 



= 2Gl\VrA''W^f3L''^- (11.5) 
The Icpton part L"'^ is 

The hadronic part is 

= Y,i-^nfS\pB - q-pxJ^{B{pB)\Ji"\X4px)) 



■{XM\4\B{pb)), (11.6) 
where Jl is the left handed current 

J2 = cr^-^-^b. (11.7) 

In the above q = Pe + Vu^ is sum of electron and anti-neutrino four- 
momentum. Wai3 is a second rank tensor. It depends on pB = uibV 
and q the momentum transfer to the hadronic system. The most general 
tensor is 

Wa0 = - gafiWi + VaV0W2 - iea0^u.V^'q''W3 

+ qaq0W4 + {Vaq0-V0qa)W5, (11.8) 

where Wj = Wj{q^ ,q-v) are scalar functions. Using this decomposition for 

+ Wsq^E, - E^Jjei-q^AEe + E^J, (11.9) 

where we have explicitly displayed the ^-function that sets the limit on 
neutrino energy. The functions W4 and W5 don't contribute to the decay 
rate since q"Lai3 = q^LajS = 0. 

The neutrino is not observed and integrating the above expression over 

E,^^ or equivalently over v ■ q = E^ + E^^ gives the measurable decay distri- 
bution dr/dq^dEg. Since q^ = 2EeE^^{l — cos9) for fixed Eg the minimum 
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value of occurs for cos 6 = 1 and the maximum for cos 9 = —1. The 
maximum neutrino energy and at fixed electron energy are, 

The maximum electron energy is at = (i.e. Ej^^ = 0), 

™ax ^ (11.11) 

2mB ^ ' 

For a given final hadronic mass mx^ , electron energy and q^ the neutrino 
energy is fixed. Using Px = {Pb — Pe — Pv^Y gives the neutrino energy as 
a function of E^, and mx^ , 

mi, + q^ — m\ 

= —5 — ^ - Ee. 11.12 

2mB ^ ' 

So integrating over E,,^ at fixed q^ and -Bg is equivalent to averaging over 
a range of final hadronic invariant masses. For values of q^,Ee near the 
upper boimdary of the {E(.,q^) Dalitz plot 

(ms - 2Ee)q^ - 2Ee{m\ - ^msEe - rrix^i^) = 0, (11.13) 

only final hadronic invariant masses near ruxmin = mo get averaged over 
in the integration over E^^ . 

The hadronic tensor W„/3 parametrizes all the strong interaction physics 
relevant for inclusive semileptonic B decay. It can be related to the dis- 
continuity of a time ordered product of currents across a cut. Consider 



^ocH = -i j d' 



4^,-i,..m4>)M^)M. (11.14) 



2m B 

Inserting a complete set of states between the time orderings 

{B\.Jljx)\X,){XAJumB) 



„ 2ms 



(11.15) 
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and performing the space-time integration gives 

At fixed q the time ordered product Tafi has cuts in the complex 
plane along the real axis. One cut is in the region — oo < < uib — 
+ rri^^rnin and other is in the region + rriL^i^ — uib < < oo. 

c y ebb 

The two cuts are well separated. The discontinuity in T across its cuts is 
evaluated using the formula 

^ P- - m6{u). 



ui + is u) 
This gives 



2m B 



(11.17) 

The first of these terms is just —Wafj- For q in the region of semileptonic 
decay the second delta function is not satisfied and so 

--ImT^p = W,,(i. (11.18) 

TT 

It is convenient to express Tq/j in terms of Lorentz scalar form factors as 
we did for W. 

+ qaqi}Ti + {vaqi3 + Vf3qa)T^. (11.19) 

The Tj's are also functions of q^ and q ■ v. One can study the T, 's in the 
complex q-v plane at fixed q^. This is the Lorentz invariant way of phrasing 
the analytic structure. For the cut associated with physical hadronic states 
containing a c quark {pB — q)—px = 0, and so m% — 2v ■ qniB + q^ = nix ■ 
Therefore 

v-q<— (11.20) 

^ 2mB ^ ^ 
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|q- V 



Fig. 2. Contour for recovering integrals of the Wj's from those of the Tj's. 



On the other hand cut corresponding to states with ebb quark content, 



{Pb + q) - Px = 0, implying that m% + 2v ■ qms + q 



V ■ q> 



{m'^x'^in -m%-q'^) 

ebb 

2mB 



, . In this case 



(11.21) 



The diflFerential decay rate dV/cPqEe is obtained from dT / dq^dEedv ■ q 
by integration over v ■ q. Integrals over v ■ q of the structure functions 
Wj{q'^, V ■ q) are related to integrals of Tj over contour shown in Fig. (2). 

The structure functions Tj can be expressed in terms of matrix elements 
of local operators using the operator product expansion (OPE) to simplify 
the time ordered product [15-18]. 



d'^xe 



'T{jUx)jL0m- 



(11.22) 



The i?- meson matrix clement of the above time ordered product is Tap. 
The coefficients of the operators that occur can be reliably computed in 
perturbative QCD in a region that is far away from the cuts. These oper- 
ators involve quark and gluon fields. They can be expanded in powers of 
l/mb by making a transition to HQET. 

Note that the contour in Fig. (2) necessarily touches the cut. However 
this is in a region that is far from the upper boundary of the (Ee, q^) Dalitz 
plot and is a place where the kinematic restrictions are on the production 
of high mass states (i.e., mx^ » ttid)- In this region the OPE is expected 
to give inclusive differential decay rates that are valid locally (i.c, without 
averaging over final state hadronic masses) up to small corrections. If the 
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contour pinched the cut near its cndpoint at the right, there would be 
very large corrections since this corresponds to a kinematic region near the 
upper boundary of the {Eg, q^) Dalitz plot where the rate is dominated by 
the production of low mass final hadronic states with masses near . 

At lowest order in perturbation theory the matrix element of this time 
ordered product between 6-quark states with residual momentum k is 

1 - (1 - 75) , , (1 - 75) 

7 2~r^'"T" — — \^bf -i + fn/3 — ^ — u. 

[mbv — q + ky — + le 2 2 

(fl.23) 

From eq. (11.23) it is evident that (without including perturbative cor- 
rections) the discontinuity across the cut in Fig. (2) is concentrated at 
its endpoint. This is another reason why it is not a problem to have the 
contour touching the cut far from this point. 



12. 6-quark Decay 

The leading contribution to the decay rate comes from the order fc" terms 
in eq. (11.23). Reducing the product of 3 gamma matrices to a sum of 
single gamma matrices it becomes 

-^u{{mbV - q)cxl0 + {mbV - q)i3^a - {mt^ - 4)9afS 

-ieafsxnirribv - q)^r}^^^Y^u, (12.1) 
where 

Ao = {rribV - qf -ml + is. (12.2) 

Matrix elements of b^i\b and tr/xj^b between 6-quark states are u'y\u and 
ujx'fsu. Hence the leading term in the OPE is obtained by replacing u 
above with 6-quark field. To get the T/s we need the B meson matrix 
elements 

{B{pB)\hxb\B{pB)) = 2pBx = 2mBVx, 

{B{pB)\hxl5b\B{pB)) = 0. (12.3) 

Note that these are exact. In this case there is no need to make a transition 
to HQET to determine the m,i, dependence of the matrix clement. The first 
of eqs. (12.3) follows from b-quark number conservation and the second 
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from parity invariance. The resulting structure functions are easily found. 
For example 

which impHes that 

TT 2 



4m6 V V 2m6 

The leading term in the OPE gives functions Wj that correspond to free 6- 

quark decay. Consequently wc have derived the usual result that at leading 
order in KQco/rnb the B meson semileptonic decay rate is equal to the h- 
quark semileptonic decay rate. 



13. The Chay— Georgi— Grinstein Theorem 

At linear order in k the 6-quark matrix element of the time ordered product 
of currents gives 

2k- {rribV-q) _ \ , r ^ 

u{{mbV - q)alp + [nibv - q)i3ja 

^0 

After converting to HQET 6-quark fields these produce terms in the op- 
erator product expansion corresponding to the operators b-u"f\iDrby and 
bvlxlbiDrby. The second of these has zero B meson matrix element by 
parity. For the first use by'yxiDrby = vxbyiDT-by. The general form for the 
B meson matrix element of this operator is 

{B{v)\byiDrby\B{v)) = Xvr. (13.2) 

Contracting with v'^ implies 

X = {B{v)\byiv ■ Dby\B{v)) = 0, (13.3) 

by the equations of motion. So there are no h.QCD/mb corrections to b- 
quark decay picture for inclusive B meson decay [16]. 
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14. Higher Order Corrections to the Inclusive B Decay Rate. 

We have seen that up to corrections suppressed by as(m(,) and Ag^^j/m^ 
the inclusive semileptonic B meson decay rate is equal to the free 6-quark 
decay rate. 

The perturbative corrections can be included by calculating the pertur- 
bative corrections to the 6-quark decay rate. This has been done to order 
alinib) [19,20]. 

At order Ag^^^/m^ dimension 5 operators occur in the OPE. They are 

the same operators that occur in the B meson mass formula. At this order 
the nonperturbative corrections to the inclusive semileptonic B decay rate 
involve Ai and A2 [17,18]. Similar results hold for the inclusive nonleptonic 
B decay rate. 

15. NRQCD 

HQET is not the appropriate effective field theory for systems with more 
than one heavy quark. In HQET the heavy quark kinetic energy is ne- 
glected. It occurs as a small l/rng correction. At short distances the static 
potential between heavy quarks is determined by one gluon exchange and 
is a Coulomb potential. For a QQ pair in a color singlet it is an attractive 
potential and the heavy quark kinetic energy is needed to stabilize a QQ 
meson. For QQ hadrons (i.e., quarkonia) the kinetic energy plays a very 
important role and it cannot be treated as a perturbation. 




Fig. 3. One loop QQ scattering in HQET. 



In fact the problem is more general than this. Consider, for example, 
trying to calculate low energy QQ scattering in the center of mass frame 
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in HQET. Setting v — Vr = (1,0) for each heavy quark ( intitial residual 
momenta are k± = (0,±k) and final residual momenta are k'^ = (0,±k')) 
the one loop Feynman diagram in Fig. (3) gives rise to a loop integral 



/ 



Q % % % % 

(27r)" (gO + ie) {-qO + ie) {{q + k+f + ie) {{q + k'^f + ie) ' ^^^'^^ 



The g° integration is ill defined because it has poles above and below the 

real axis at ~ zLie. This problem is cured by not treating the heavy quark 
kinetic energy as a perturbation but rather including it in the leading order 
terms. 

Properties of quarkonia arc usually predicted as a power series in v/c 
where v is the magnitude of the relative QQ velocity and c is the speed of 
light. For these systems the appropriate limit of QCD to examine is the 
c — » oo limit, where the QCD Lagrangian becomes an effective field theory 
called non-relativistic quantum chromo dynamics (NRQCD) [21]. For finite 
c there are corrections suppressed by powers of 1 /c. In particle physics we 
usually set Planck's constant divided by 27r and the speed of light to unity. 
Making the factors of c explicit the QCD Lagrangian density is 

Cqcd = - Jg^.G^"^ - cQ{jp - mQc)Q. (15.2) 
In the above the component of a partial derivative is 

9o = \l (15.3) 
and D is the covariant derivative 

D^ = df, + ^A^T^. (15.4) 

The gluon field strength tensor G^^ is defined in the usual way except that 
9 g/c- 

All dimensionful quantities can be expressed in units of length [x\ and 
time [t]. (i.e., [E\ ~ l/[t\ and \p] - l/[x\). The Lagrangian L = j dPxC 
has units of l/\t] since the action S = J Cdt is dimensionless. It is straight- 
forward to deduce that the gluon field has units [A\ <^ l/^J[x\\t] and the 

strong coupling g ^ ^ [x] / [t] . The fermion field has units [V-j] ^ l/[.x]^/^ 
while its mass has units [toq] ~ M/N^- With these units mqi? has di- 
mensions of energy and the strong fine structure constant a = /A-kc is 
dimensionless. 

For the fermion field Q the transition from QCD to NRQCD follows the 
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derivation of HQET. It is rewritten as 



1 + 

WLQC 



(15.5) 



where ^ is a two-component Pauli spinor and the _L is with respect to 
the rest frame four- velocity Vr = (1,0). Using this, the part of the QCD 
Lagrange density involving Q becomes 

72 



+ . . . , (15.6) 



where the ellipses denote terms suppressed by powers of 1/c. Note that the 
heavy quark kinetic energy is now leading. The replacement g ^ g/c was 
necessary to have a sensible c — > cx3 limit. 

Among the terms suppressed by a single power of 1/c is the gauge com- 
pletion of the kinetic energy 

Ci„t = ^A^[V+T^V^ - (V^)tT^V]- (15.7) 

mqc 

There is also a 1/c term involving the color magnetic field = V x A. 

It is convenient to work in Coulomb gauge, V • = 0. Then the part 
of action that involves the gluon field strength tensor and is quadratic in 
the fields simplifies to 

= ^J (fx {{diAj^f + {doAff - {diAff) + nonabelian terms. 

(15.8) 

The component of the gauge field does not propagate, and its effects on 
the heavy quarks is reproduced by an instantaneous Coulomb potential. 

K=(x,y) = |^^^. (15.9) 
4tt |x — y| 

In QQ systems the component of the gluon field typically has a four- 
momentum of order {E,p) ~ {rnqv^ ,mQv). It is far off-shell and can 
be integrated out of the theory giving rise to the Coulomb potential 
above. The transverse gluons A have both potential modes where {E, p) ~ 
{iTiQV^ ,mQv) and propagating modes, where (£^, p) ~ {mQV^ , rnqv^ / c) , 
that are important [22]. The potential modes do not propagate and as 
was done for Aq the Lagrangian can be written without these fields ex- 
plicitly appearing. However the propagating modes must appear explicitly. 
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For these modes the gluon field is rewritten in terms a scaled spatial com- 
ponent y = x/c so that the {1 / 0^)8^ / dt"^ part of the kinetic term is as 
important as the part of the kinetic term [23] , 

A^(x,i) = ^A^(x/c,f). (15.10) 

Using this the QCD Lagrangian becomes, for large c, 

Lqcd = Lnuqcd + ■ ■ ■ , (15.11) 
where the ellipsis denote terms suppressed by powers of 1/c and 

- j rf^a^id^xaV Vi,t)r^V(xi, i)K(xi, X2)V^ (x2, t)r^V(x2, i) 

4/''°'((l-*?(>'-'))'-(i|^-^?(^-'))')' f''^-'^^) 

with y — x/c. 

At leading order in the 1/c expansion the radiation gluons do not in- 
teract with the heavy quarks. Also the nonabelian terms in G^^G^^" are 
suppressed. The leading interaction of the propagating gluons with the 
heavy quark is 

^ir^t = ^^^^A^(0,i)[VtT^VV - (V^)tT^V]- (15.13) 

Note that the term involving the propagating gluon color magnetic field is 
less important because Be = V x A (x/c, i) = + ^ V x A(0, t) + . . . is 
suppressed by 1/c. This derivation of NRQCD has followed closely that in 
Ref. [23]. (For a somewhat different scaling of fields see [24].) 

The NRQCD Lagrangian does not have a heavy quark flavor symmetry 
because niQ explicitly appears in the kinetic energy term. However it 
does have the spin symmetry. Furthermore, the leading interaction of the 
propagating gluons with the heavy quarks also respects the heavy quark 
spin symmetry. 

In some respects the Feynman rules for this theory are a little unusual. 

Coulomb exchange is denoted by a dashed line corresponding to a propa- 
gator — z/k^ and a vertex —igT^. The heavy quark, tp, propagator is 
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The transverse radiation gluons are denoted by the usual wavy hne. Their 

propagator is 



(15.15) 



(fc0)2-k2+i£' 

and their interaction with heavy quarks is determined by eq. (15.13). It 
gives rise to Feynman rules where energy is conserved at the vertices but 
not the three-momentum. The transverse radiation gluons do not transfer 
any three -momentum to the heavy quarks. 

Anti-quarks can be added in the same way as in HQET. For these one 
writes 



1 + 



2mQC 



+ . 



(15.16) 



and their coupling is similar to that of the quarks. The field x creates a 
heavy anti-quark while destroys a heavy quark. 



Fig. 4. Contribution to matrix element in Eq. (15.17). 



The formalism developed above helps us understand the size of matrix 
elements in the v/c expansion. For the first example consider 

(35i|(VVx)(xVV)|'5i). (15.17) 

Here | ^Si) denotes a QQ state at rest with quantum numbers (with the 
usual (^•S'+i)^; spectroscopic notation). For cc quarkonia the lowest mass 
state with these quantum numbers is the if) and for bb it is the T. At leading 
order it is only the Coulomb potential gluons that are relevant. Feynman 
diagrams where a potential glnon is exchanged between the quark-anti- 
quark pairs in eq. (15.17) vanish. Fig. (4) is such a diagram. The loop 
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integration is of the form, 

1 1 i 



(15.18) 

which is zero because all the poles in the fc" integration are on the same 
side of the real axis. This generaUzes to all orders in perturbation theory 
and consequently the matrix element factorizes. 

( 35i|(vV^x)(xV'V)l 'Si) = ( 35i|(vV^x)|o)(o|(xVV)| '5i). 

(15.19) 

Corrections to factorization come from radiation gluons and so are sup- 
pressed by powers oiv/c. 

Next consider a color-octet matrix element, 

( 35i|(V'^T^a'x)(xV^r^V)l '5i). (15.20) 

Since | Si) is a color singlet at least one radiation gluon must be exchanged 
between the quark anti-quark pairs. Such a diagram is shown in Fig. (5). 
But the coupling of the transverse gluon is proportional to p and averaging 



Fig. 5. Contribution to matrix element in Eq. (15.20). 

with the ^-wave, | ^S\) momentum space wavefunction, causes this diagram 
to vanish. Two transverse gluon exchanges are needed to get a non-zero 
ampUtudc and so this matrix element occurs at order (y/cf'. (Recall £j„t 
is proportional to 1/c^/^.) 

This power counting is correct for very heavy quarks where rnqv^ /c ^ 
^QCD- For bb and cc systems this inequality is not satisfied. Often for 
these systems one uses a modified power counting which is appropriate for 
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rnqv^ /c ^ ^qcd- In this situation the coupling of the radiation gluons 
should be nonperturbative with a of order unity and g ~ y/c. Then Lint 
is suppressed by only 1/c and the leading contribution to the color-octet 
matrix clement in cq. (15.20) is suppressed by (w/c)** instead of {v / cf" . 

It is easy to see that with ~ -^c for radiation gluons their nonabelian 
self couplings are not suppressed. For example, the three gluon coupling 

tC-igiuon ~ d^x {VA(x/c,t))A(x/c,i)A(x/c,i)} 

~ ^ y d?y{V yA{y, t))k{y, t)A{y, t), (15.21) 

is of order unity if 5 ~ \/c. Since the quark coupling to the radiation gluons 

is suppressed, even with raQvjc? ^ Kqcd, the leading order quarkonia 
states are still determined by the Coulomb potential. (Note that the strong 
coupling for potential gluons retains the usual scaling 5 1.) For the rest of 
this section this modified power counting appropriate to mgw^/c ^ Aqcd 
is used. (It would be useful to have a more systematic understanding of 
this modification of the power counting.) 
Next consider the matrix element 

( 35,|(^t_^TVx)(x^^rVV)i'5i). (15.22) 
mgc ruQC 

The factors of mgc in the denominator arc inserted so that this operator 
has the same dimension as those occuring in eqs. (15.20) and (15.17). Now 
only a single radiation gluon is needed between the quark pairs. This matrix 
element is also suppressed compared to that in cq. (15.17) by (w/c)^. Two 
of the factors of 1/c coming from the radiation gluon vertex and two from 
the explicit c's in the operator. 

Radiation gluons have a novel coupling to heavy quarks. Usually in the 
derivation of the Fcynman rules there is an integral J d'^k/(2Tr)^ for each 
propagator and a {2TT)^S^{ki — kf) for each vertex. The result is that there 
is one / (i''fc/(27r)^ integration for each loop. But the leading vertex for the 
ip[p) 'ip{p')A{k) radiation gluon heavy quark interaction in Lint contains 
the delta function {2'k)^6{E - E' - k'^)5^{p - p'). At leading order in v/c 
the radiation gluons do not transfer three-momentum. One consequence of 
this is that the matrix element 

(35i|(Vt-^TVx)^...^(xt-^rVV)| (15.23) 
ruQC ruQC mqc^ mqc 
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is suppressed compared with 

( 35i|(V>t-^TVx)(x^^^ • • • ^tVV)| 'S,). (15.24) 

mgc mQcniQC mqc 

For the former matrix element the Feynman diagram in Fig. (5) contains 
the factor {p' — pY^ . . .{p' — p)'"" from the spatial derivatives on the QQ 
quark pair. In the rest frame of the "^Si state this vanishes unless one 
expands A{x/c,t) in £j„t to m powers of 1/c [25]. The m'th term in 
the expansion of A gives a factor x'^^ . . . x''"/c™9ii . . . di^A{0, t) producing 
in the Feynman diagram the factor (fc*^ . . . k^^) / di^ . . . di^5'^{\)' — p). 
(The factors of Xt become derivatives di acting on the delta function.) The 
Feynman integration over p' is done by integrating by parts m times 

/ dh'{p' - p)j, ...{p'- p),^di, . . . dij\^' - p) 



= (5iiji . . . 5i^j^ + perms). (15.25) 

The matrix element with the m derivatives acting on the pair is sup- 
pressed by (w/c)™ compared with the matrix element where the m deriva- 
tives act between the quarks in the x^V' pair. 

One of the most important recent developments in quarkonia physics is 
the realization that octet matrix elements suppressed by powers of v jc are 
sometimes more important than the analogous color singlet matrix element 
because they have coefficients that are enhanced by factors of l/as(mg). 
A simple example of this occurs for production of quarkonia (with 
quantum numbers) at large momentum transverse to the beam in pp scat- 
tering [26]. At the present time this has only been observed for the case 
Q = c corresponding to V or tf)' production. 



Fig. 6. Color-octet quark anti-quark production from fragmenting gluon. 

The production of quarkonia at large transverse momentum is dominated 
by gluon fragmentation. The underlying process is a virtual gluon k"^ ~ 
4to^ converting into a cc pair in a '^Si, color octet configuration. The 
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color octet cc pair omits radiation gluons producing a "^X final state. The 
differential cross sections is written as 

da = J dzdagg^gg{kg = p^/z)Dg^^{z), (15.26) 

where the ^ carries a fraction z of the gluons momentum. The square 
of the diagram in Fig. (6) gives a gluon fragmentation function that is 
proportional to 6{1 — z). Explicitly 

Dg^^{z)<xg\2m,)(^-^^ 5(1 - ^)(0('^^)), (15.27) 

where (0*^ '^^^) is proportianal to the fragmentation probability. In pertur- 
bation theory it can be written as, 

(O(^^^i)) =^(0|(xlTVV'c)|^ + ^)-(^ + ^|(^]rVxc)|0). (15.28) 

X 

The hadronic matrix element above is evaluated in the '4> rest frame and 
a sum over i]j polarizations is understood. The fragmentation function is 
proportional to 5{1 — z) since the radiation gluons carry away a fraction 
of the ■!/''s momentum suppressed by (w/c)^. The factor of (1/4to^)^ in 
eq. (15.27) comes from squaring the gluon propagator. To write the cross 
section in the form above, the phase space factor d^p^/2p''^ is set equal to 
d^ kg /2kg since for a high energy process the ^ mass can be neglected in 
p^. The radiation gluons hadronize into a set of pions accompanying the 
i/j that carry a fraction of order iP' j c? of the V-''s energy. 

The fragmenting gluon is almost on shell. An on shell gluon has only 
transverse polarizations and so the color octet cc pair produced by the 
fragmenting gluon is also transversely aligned. Since the leading couplings 
of the radiation gluons to the heavy quark preserve spin symmetry this 
transverse alignment is transferred to the in the final state [27]. (There 
are vjc [28] and perturbative ols corrections [29] that reduce this alignment 
somewhat.) 

There is a color singlet contribution to the g ^ ip fragmentation func- 
tion. It is calculated from Feynman diagrams like Fig (6), but with two 
hard gluons radiated off the quark legs to make the heavy quark anti-quark 
pair in a color singlet configuration. This contribution to the fragmentation 
function is enhanced to (c/w)^ compared to that in eq. (15.27). However, it 
is suppressed by (as(2mc)/7r)^ because of the two additional hard gluons. 
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16. Conclusions 

Over the last decade there has been remarkable progress in our ability to 
predict properties of hadrons containing a single heavy quark and proper- 
ties of hadrons containing a heavy quark anti-quark pair (i.e., quarkonia). 

The application of heavy quark symmetry and the operator product ex- 
pansion allows for model independent predictions for exclusive and inclusive 
B decays. These predictions will play an important role in the determina- 
tion of the 6 — > c and b ^ u elements of the Cabibbo-Kobayashi-Maskawa 
matrix. These lectures provide a rudimentary introduction to these meth- 
ods. 

Chiral perturbation theory and Lattice Monte Carlo techniques are also 
useful in conjunction with the heavy quark methods developed in my lec- 
tures. The latter has been discussed in other lectures in this school and 
reviews on combining heavy quark and chiral symmetries can be found in 
Refs. [30,31]. 

Recent advances in NRQCD have improved our understanding of the 
decays and production of quarkonia. These results are particulary dramatic 
for the the production of quarkonia at high energy accelerators where the 
octet matrix elements often dominate. 

NRQCD and HQET remain active areas of study and development. 
Even without the discovery of dramatic new theoretical techniques the 
continuing experimental developments will provide new and challenging 
opportunities to apply these effective field theories. 
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